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Abstract. Crystalline silicon is currently being discussed as test-mass material 
for future generations of gravitational wave detectors that will operate at cryogenic 
temperatures. We present optical absorption measurements on a large-dimension 
sample of crystalline silicon at a wavelength of 1550 nm at room temperature. The 
absorption was measured in a monolithic cavity setup using the photo-thermal 
self-phase modulation technique. The result for the absorption coefficient of 
this float-zone sample with a specific resistivity of llkficm was measured to be 
OA = (264 ± 39) ppm/cm. 



1. Introduction 

The initial (1st) and advanced (2nd) generations of interferometric gravitational wave 
detectors employ suspended fused silica test masses and use a laser wavelength of 
1064 nm [T]-[3]. Observatories beyond the 2nd generation will require very high laser 
powers to reduce the quantum noise at frequencies above ^ 50 Hz, while at lower 
frequencies it is promising to cool the test masses to cryogenic temperatures to reduce 
thermal noise. 

While at room temperature fused silica shows a high mechanical Q-factor [3] 
and low optical absorption [5], the mechanical Q-factor decreases by several orders 
of magnitude at cryogenic temperatures [SlIZ]- This makes fused silica unsuitable as 
test mass material for cryogenically operated GW observatories. Crystalline silicon, 
however, shows a promising mechanical Q-factor at room temperature that even 
increases towards cryogenic temperatures up to 2 x 10^ [HIS]- Due to the very high 
absorption coefBcient of about 10 /cm at 1064 nm, silicon test- masses require a change 
to higher laser wavelengths, where the absorption coefficient decreases rapidly jlOj . 

A wavelength of 1550 nm is located within the silicon energy gap and appears 
promising because of the availability of appropriate optical components and lasers 
developed for the telecommunication sector. Also, strong squcczcd-light sources of 
> 12 dB arc available at 1550 nm [11] to increase the sensitivity of GWDs beyond 
the quantum limit |12j . An up to now open question are measurements of the optical 
absorption coefficient at 1550 nm in the temperature range from a few kelvin up to 300 
K. Unfortunately, the measurements presented in |10j were done at slightly shorter 
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wavelength up to 1450 nm. Furthermore, these measurements are a prediction of the 
absorption coefficient that is based on photo-current measurements. They might not 
include various effects that lead to optical absorption but do not generate charge 
carriers. Also the photo-current measurement might not show all absorption effects 
that are relevant for thermal noise. 

In this paper we present direct optical absorption measurements on two silicon 
samples in large dimensions at room temperature using the photo-thermal self- 
phase modulation technique |13j . Silicon sample A forms a monolithic cavity. The 
measurement technique is perfectly adapted to this setup and therefore gives very 
precise results. While the photo-thermal effect delivers the absorption coefficient, 
with this method the round-trip loss is measured independently at the same time. 
Since losses apart from absorption are small for this monolithic setup, the two results 
confirm each other. In a second experiment, an uncoated silicon sample B from 
another manufacturer was investigated to ensure that sample A was not coincidentally 
untypically high or low absorbing or that a large amount of the absorption was not 
caused by the coating process or by the coating itself. The uncoated sample was placed 
in a Fabry-Perot cavity. This measurement only served as an order of magnitude 
estimation that yielded an approximate upper and lower limit of the absorption and 
entirely confirmed the result of the first measurement. The slightly different setup 
additionally backed the fact that the optical absorption measured for sample A was 
not mainly caused by the dielectric coatings. 

2. Absorption Measurements on Silicon Sample A 

Sample A was manufactured by Siltronic AG |14j with the Czochralski technique. 
The crystal's orientation is (111). The material has a specific resistivity of 
about llkficm, which indicates a low doping or contamination with foreign atoms. 
According to the manufacturer it is a low boron doping, which is a p-donator. This 
means the impurity concentration is approximately 2 x 10^^ atoms per cm"^ |15j . 
Sample A is the purest material available to us in the required dimensions at ordering 
time. 

The substrate was cut and polished into a cylinder with the rotation axis being 
parallel to the (111) axis. The cylinder's diameter was 2 x R = 10 cm, the length 
was i = 6.5 cm. The end surfaces were polished to be convex curved with a radius of 
curvature of 1 m to form a cavity with a free spectral range (FSR) of 663 MHz. 

The substrate's curved end surfaces were coated using ion beam sputtering (IBS). 
The high-reflection coatings consisted of Si02 and Ta205 and had a design reflectivity 
of > 99.9% at a wavelength of 1550 nm. Hence, the coated substrate formed a 
monolithic cavity with beam propagation along the (111) axis. 

2.1. Experimental Setup 

Fig. [T] shows a schematic of the experimental setup. A laser beam at a wavelength of 
1550 nm was mode-matched to the cigenmodc of the monolithic cavity. To calibrate 
the time-axis of our measurements we used frequency markers. An electro optical 
modulator (EOM) generated these frequency markers by imprinting sidebands at a 
frequency of 43.57 MHz onto the light field. 

Photo detector PDi detected the reflected light, which was separated from the 
incident field by a combination of a Faraday rotator and a polarizing beam splitter 
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Figure 1: Schematic of the experimental layout. The cylindric silicon substrate 
with convex curved and coated end surfaces formed a monolithic cavity. A function 
generator (FG) that actuated the piczo electric transducer (PZT) of the laser 
modulated the laser wavelength. Photo detector PDi detected the power reflected 
by the cavity and showed the cavity resonance peaks. On the right side, a photograph 
of the monolithic cavity with 6.5 cm length and 10 cm diameter is shown. 

(PBS). Demodulating PDi and creating a Pound-Drever-Hall-type [16] error-signal 
generated the frequency markers. 

For exploiting the photo-thermal self-phase modulation technique, the laser 
frequency was scanned around the resonance frequency of the cavity via a piezo electric 
transducer (PZT). An increasing wavelength corresponded to a shortening of the cavity 
and a decreasing wavelength to a lengthening of the cavity. 

The input laser power was 22 mW for all measurements. We performed several 
measurements by varying the scan frequency starting from 0.2 Hz to 49 Hz in 18 
steps with a constant scan amplitude. Because the PZT showed a hysteresis, the 
actual wavelength change had to be calibrated for each frequency for increasing and 
decreasing wavelength. The time axis of each measurement was calibrated from the 
scan frequency, the FSR, and the detected error signal. In Figure [2l an example 
measurement for a scan frequency of 0.6 Hz is shown. This frequency corresponds to a 
scan velocity of about 2ms/peak at half width half maximum (HWHM). The yellow 
crosses (narrow peak) show the measured peak for a decreasing wavelength with the 
corresponding simulation in red (solid line). The light blue crosses (broad peak) show 
the peak for an increasing wavelength with the simulation in dark blue (solid line). 
Without absorption, the two peaks were identical. 

2.2. Measurement Analysis and Results 

To calculate the absorption coefficient a from the measured peaks as shown in Fig. 
the peaks were Jittcd. For the fitting process we used the parameters from Table [T] 
as well as Ri, i?2;^nd for measurements with a visible thermal effect a as fitting 
parameters. Here, i?2 is the effective reflection of the out-coupling coating accounting 
for the entire cavity round trip loss apart from the transmission of the in-coupling 
coating. The values for material parameters were taken from literature |17fll9| . For 
the geometric parameters, we used values based on our best knowledge of the cavity 
design. A Nelder-Mead algorithm was run to find the best set of fitting parameters 
minimizing the standard deviation of measurement and simulation. 13 of 18 single 
measurements showed a visible thermal effect and were used to derive the absorption. 
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Figure 2: Example of measured (crosses) and simulated (solid lines) reflection peaks at 
a scan frequency of 0.6 Hz: The broad (blue) peak forms for an increasing wavelength, 
the narrow (orange) one for a decreasing wavelength. Without absorption, both peaks 
would be identical. 



Table 1: Material and geometric parameters of the monolithic silicon cavity used for 
the simulations. 



Geometric parameters 




Source 


crystal radius R 


5 cm 


specified by manufacturer 


crystal length L 


6.5 cm 


specified by manufacturer 


beam waist 


217.8 ^m 


calculated from L, n and ROC 


Material parameters 


index of refraction n 


3.48 


m 


thermal expansion oth 


2.62 X 10- VK 


m 


thermal refr. coefF. dn/dT 


1.87 X 10"'^/ K 


m 


specific heat c 


714J/(kgK) 


m 


density p 


2330kg/m3 


m 


thermal conductivity k^h 


140 W/(mK) 


m 


material emissivity e 


1.0'^ 





"0.0 < e < 1.0 are the boundaries for the thermal emissivity. For our systems the 
value of this parameter is not relevant since R^ luq- 



The remaining five measurements were used to derive Ri and i?2 and served as a 
consistency check. 

Figure [DJa) shows the results for the absorption coefficient a derived from the 
13 different measurements (dark-green dots). The purple lines show the mean value 
of all single results and their standard deviation which is a = (264 ± 39) ppm/cm 
(39ppm ^ 15%). 

The results for the (power) reflection Ri of the in-coupling coating are shown in 
Figure [3ljb) . The mean value and standard deviation of all 18 single measurements 
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Figure 3: (a) Results for the absorption from single measurements (green dots). The 
mean value and the standard deviation of a = (264 ± 39) ppm/cm are given by the 
purple line and the dashed purple lines, respectively. The dots in (b) and (c) show the 
results for the in-coupling reflectivity Ri and the effective out-coupling reflectivity R2 
with Ri = (99.9784 ± 0.0015) % in orange and R2 = (99.630 ± 0.025) % in turquoise. 
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is i?i = (99.9784 ± 0.0015) %. Measurement numbers 14-18 did, as stated above, not 
exhibit a visible thermal effect. 

Figure |3Jc) shows the results for i?2 for all single measurements (dark-green 
dots) and their mean value and standard deviation (turquoise lines) of i?2 = 
(99.630 ± 0.025) %. This result was used to cross-check the obtained absorption 
values. An absorption coefficient of a = (264 ± 39) ppm/cm results in an absorption 
of atotai = 264 ppm/cm x 13 cm = 0.343% per round-trip. The round-trip loss added 

to i?2 results in a new effective reflection i?2 = i?2 + Q^totai — 99.973% that still 
contains the cavity scattering loss. Since the pure reflection can be assumed to be 
very similar to i?i (identical coating design, but different coating runs) the results for 
the three parameters agree perfectly. This is an additional consistency check and not 
an automatical consequence of the simulation. 

We repeated the series of measurements for a polarization rotated by 90° . Further 
we exchanged in-coupling and out-coupling coating and repeated the measurements 
for the two polarizations. As expected, in each case the results for the absorption 
coefficient agreed with the result presented above within the error bars. 

2.3. Error Propagation 

The error bar of 39ppm or 15% corresponds to the standard deviation of 13 
independent measurements using different scan-velocities. An additional error bar 
arises from uncertainties in the simulation input parameters. To estimate this error, we 
individually changed the input parameters listed in Tabled] by ±10 % and recalculated 
Ri, i?2 and a. We found that for most parameters the influence on a is approximately 
linear and none of the changed parameters caused a change of the result for a by more 
than 15 %. 

For n, the error bar of the value from literature is in the order of 10~^ |17j and 
therefore negligible, the uncertainty of dn/dT is in the same order of magnitude |19j . 
ath and dn/dT affect the result as a sum. Since oth ^ dn/dT, an uncertainty of Oth 
is negligible. The error bars of c, p and fcth are not known to us. We estimate that the 
uncertainties of the cavity geometric input parameters as well as for the calibration of 
the time axis and in the measurement of the mode-matching are below 10%. 

If the simulation input parameters are precise within 10%, our statistical error 
bar of 15% is the dominating error contribution. 

3. Measurements on Silicon Sample B 

The absorption measured with sample A was unexpectedly high (see discussion). To 
verify that this sample did not absorb untypically much, or that the coating process 
or the coating itself did not cause the absorption, a second experiment with a different 
sample B was performed. The two samples differed in manufacturer and crystal 
orientation. 

A cylindric substrate with one inch diameter and 3.5 cm length was manufactured 
by Mateck [20j in (100) orientation using the Czochralski procedure. The (100) axis 
is the rotation axis of the cylinder. The cylinder end surfaces were polished to be 
parallel with a tolerance of < 30". The material was declared by the manufacturer to 
be undoped with a resistivity of > 5 kilcm. 

The substrate was placed in a resonator at Brewster's angle of 74 ° to minimize 
round trip loss due to reflection at the substrate surface (see Figure |4]). The mirrors 
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Figure 4: Experimental setup for the measurement on the second silicon sample. The 
beam was coupled into the cavity through mirror Mi. The silicon substrate was placed 
in the cavity at Brewster's angle of 74 ° which caused a strong beam displacement. M2 
was the end mirror of the cavity. The resonance peaks were detected in reflection of 
the beam splitter (BS). The laser wavelength was scanned with a function generator. 

were clamped to an aluminium spacer. The laser beam coupled into the cavity through 
the in-coupling mirror Mi. The reflected beam power Prcfl was detected in reflection 
of the beam splitter (BS). The modulation of the laser, the calibration of the time 
axis and the measurement procedure were identical to the procedures described in 
Section UA\ for the monolithic cavity setup. 

Entering the substrate at Brewster's angle causes an elliptical beam profile within 
the substrate. Since not necessary for other experiments, a non-circular beam profile is 
not implemented in our simulation program. The discussed measurement had the goal 
of independently determining a lower limit for the optical absorption to confront the 
obtained value with the results from sample A. Therefore, a calculation of an upper 
limit for the power density and thus for the heat distribution within the substrate 
is sufficient. This is given by a circular beam profile with the radius of the minor 
semi-axis of the elliptical profile. Since a higher power density requires a smaller 
absorption to cause the same thermal effect, this assumption yields a lower limit for 
the absorption coefficient. 

From 52 single measurements, the lower limit for the absorption was found to be 
as — (149±79) ppm/cm. The results of the single measurements arc shown in Figurc[5] 
(green dots). The mean value of all measurements with the standard deviation are 
depicted by the light-blue line and the light-blue dashed lines, respectively. Despite 
a shared spacer for mirrors and substrate, the external cavity setup proved to be 
instable and prone to acoustical disturbances, which caused large error bars due to 
the statistical fluctuations of the detected peaks. 

The error propagation was already discussed for sample A. The uncertainty of the 
result is dominated by the large standard deviation, while the errors in the material 
parameters are negligible in first order approximation. 

Apart from the elliptical beam profile, the large statistical error shows that 
the second setup is much more instable and therefore disadvantageous compared 
to the monolithic setup. Nevertheless, the lower limit for the absorption of ae = 
(149 ± 79) ppm/cm allows the conclusion that the absorption in sample A did not 
primarily originate in the dielectric coatings. 

In a second series of measurements, the round trip loss of the cavity was measured. 
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Figure 5: The results for the lower limit of the absorption coefficient of sample B 
obtained using photo-thermal self-phase modulation are shown by the green dots. 
The light blue line illustrates the mean value of all single measurements of ae — 
149ppm/cm and the dashed light-blue line the error bars of 79ppm/cm. 

This provided an the upper limit for the absorption coefficient. For this measurement, 
two mirrors with identical coatings were used and the round-trip loss was minimized 
by inclining and rotating the substrate. The maximum impedance matching was found 
to be (21.2±0.3) %. (An impedance matching of 100 % means that the reflected power 
at resonance is zero.) Using the design reflectivity of Ri=:R2 ~ (99.97 ± 0.01) % the 
optical loss was calculated to Lrt = (4400±1200) ppm. Lrt contains the entire optical 
loss that consists of absorption of and scattering at the mirror coatings as well as of 
the reflection and scattering Pout at the silicon substrate surface. The latter occurs 
twice per round trip because of entering and leaving the substrate, respectively. These 
reflections are caused by non-perfect plan-parallel end surfaces of the substrate, the 
wave-front distortion of the beam, scattering and limitations in fine-adjustment. The 
laser beam passes 2 x 3.5 cm/cos(90° — 74°) = 7.28 cm of the substrate per round trip. 
This results in a loss of (4400 ± 1200) ppm/7.28cm (604 ± 165) ppm/cm and forms 
the upper limit for the absorption coefficient. 

Nevertheless, the result of 70 ppm/cm < a < 770 ppm/cm for silicon sample B 
suggests that the result obtained with sample A was typical for samples of the degree 
of purity involved. 

4. Discussion 

Fig. [S] shows the absorption result from the present work in comparison to earlier 
absorption results from spectral response measurements on solar cells taken from 
literature (purple dots and green triangles) [2T1[22] and to a theoretical prediction 
of the absorption caused by the residual boron contamination of sample A [23j . 

Previous measurements |2HI22| found a much lower absorption of silicon than 
measured here. However, these results are also not consistent. Although no 
measurements at 1550 nm are available, the absorption measurements from Keevers 
and Green predict an absorption coefficient of a < 0.02 ppm/cm at 1550 nm |10j . while 
the measurements from Anagnostopoulos predict an absorption coefficient in the order 
of a ~ 50 ppm/cm 22\ . Green and Keevers explain their much lower absorption results 
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Figure 6: The purple dots and green triangles show band-band absorption values Q!bb 
of silieon from literature [T0j[2T] (purple dots) and [22] (green triangles). The free 
carrier absorption apc for a p-doping of iV = 2 x 10^^ /cm'^ (sample A) was calculated 
using [23] (blue line). The red dot shows our result for sample A. 



by suggesting contamination of Anagnostopoulos' sample, due to doping or unintended 
foreign atoms. The assumed kind of contamination is, however, not specified in either 
of the two publications. The band-band absorption asB is the lower absorption limit 
for intrinsic crystalline silicon. 

Even very low doping or contamination of the silicon samples can dominate the 
band-band absorption in the infrared region due to the free carrier absorption apc |23j . 
Since the semiconductor industry is by far the largest area of application for crystalline 
silicon, samples almost always are doped or at least slightly contaminated with the 
doping material that is normally used in the apparatus for the crystal growth. A 
specific resistivity of 11 ki7cm, which is the specific resistivity of sample A, corresponds 
to a p-doping of TV = 2 x 10^-^ /cm"^ (N is the number of doping atoms). The number 
of doping atoms was calculated from the specific resistivity using [15] . The blue line in 
fig.[6]shows apc for a p-doping of = 2 x 10^^ /cm'^ (blue line) calculated using [23] . 
apc Eirid the number of doping atoms are proportional |24| . This theory predicts an 
absorption coefficient of a « 1 ppm/cm for sample A, which is two orders of magnitude 
below our measurement results. Therefore, this doping theory does not explain the 
absorption of samples A and B. However, the theories for absorption due to doping 
generally are optimized for much higher doping than in our case. It is therefore 
possible that the prediction for residual-doping-induced absorption is not accurate, 
and further theoretical and experimental investigations are required to clarify this 
issue. Another possible and plausible explanation for the deviation from the numbers 
published by Keevers and Green is that their photo-current measurements did not 
include effects that, while leading to optical absorption, would not lead to generation 
of carrier charges. Such photo-current measurements are thus not unconditionally 
transferable to optical absorption and may rather be used to derive a lower limit. 
Finally, the specific resistivity values provided by the manufacturers not necessarily 
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provide information about all kinds of contamination present in the crystal and thus 
can not be used as an absolute measure for the crystal's purity. Further investigation 
of the residual contamination is required to derive a model for the dependence of the 
optical absorption on residual doping and contaminations, which are not reflected in 
the specific resistivity value. 

5. Conclusion 

We measured the room temperature absorption coefficient of two silicon samples 
A and B which differed in manufacturer and crystal orientation. According to 
the manufacturer, sample A had a residual contamination by boron atoms with 
N = 2x 10^^ /cm'^. For sample B, the specific resistivity was specified to be < 5 kficm, 
this corresponds to a p-doping of < 4 x 10^^ /cnr^. 

For sample A, a result of q;a = (264 ± 39) ppm/cm was obtained. The consistent 
results for optical absorption and round trip loss induce that no other process apart 
from the optical absorption provides a significant loss contribution. 

For sample B, upper and lower limits of 70 ppm/cm < ae < 770 ppm/cm were 
derived for the absorption coefficient. This confirmed the result for sample A not 
to be extraordinarily high or low. Since uncoatcd, the measurement of sample B 
additionally proved that the absorption of sample A did not originate in the dielectric 
coatings but was due to absorption in the bulk substrate. 
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